Low-mass quark stars as Maclaurin spheroids 
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Abstract 

We show that in the limit of low masses (M < O.IM©), the orbital and 
rotational frequencies in exact relativistic numerical models of rotating, ax- 
ially symmetric, quark stars coincide with those for Maclaurin spheroids. 
In particular, when the eccentricity of the spheroid exceeds a critical value 
(e > 0.83458318), circular orbits in the equatorial plane are unstable for a 
range of orbital radii outside the stellar surface. The orbital frequency in 
the marginally stable orbit, and all other orbits, around Maclaurin spheroids 
goes to zero in the limit e = 1. Maximum orbital frequencies in accreting 
Galactic X-ray sources may be a poor indication of the source mass. 
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As noticed by Thomas Simpson in 1743 and by d'Alembert, the rotational 
frequency of a Maclaurin spheroid (Fig. 1) 

= 27rGp{l - e2)V2e-3 {(3 - 2e^) arcsine - 3e(l - e^)^/^} (1) 

approaches zero not only for e — > 0, a "spheroid which departs only slightly 
from a sphere," but also for e ^ 1, a "highly flattened spheroid" [1]. 

We find [2] that the latter property is shared by orbital frequencies, u!{r). 
Specifically, at the equator, 

cu\a) = 2TrGp{l - e^/^e-^ {arcsin e - e(l - e^f/^} . (2) 

In fact, the corresponding circular orbits (of radius r = a) are stable only 
for eccentrities e < = 0.83458318, but the frequency of the innermost 
(marginally) stable circular orbit, present in the equatorial plane at r > a 
for all e > 65, also goes to zero for e — > 1, 

o;^, = 0.52 76189 x 27iGp{l - e^fl'^e'^. (3) 

The expressions (2), (3), for the maximum frequency in circular orbits are 
plotted in Fig. 2. Here, e = (1 — is the eccentricity of the spheroid, 

with a the major and 6 the minor axis of the same. The radius of the 
marginally stable orbit is r = r^s = 1.19820292ae > a. 

The above expressions are valid in Newtonian gravity, where the Maclau- 
rin spheroids are bodies of uniform density in hydrostatic equilibrium. In 
general relativity there arc no incompressible bodies, but quark matter is 
very nearly so, with the speed of sound close to c/3, one third the speed of 
light. 

We have computed numerical models of rotating quark stars [3] in gen- 
eral relativity and have found that for stellar masses much less than that 

of the sun (M << Mq), the density is nearly uniform throughout the star 
and, as expected, close to the density po at zero pressure of quark matter 
(assumed to be self-bound). The stars were constrained to be axisymmetric. 
We have found that for M < O.IMq the Maclaurin spheroid is a very good 
approximation to the quark star, at least in terms of gross stellar properties 
(Figs. 1, 2). 

In Fig. 3 we show the results of our general-relativistic numerical compu- 
tations of the maximum orbital frequencies, / = (27r)~^a;, in stable circular 



2 




Figure 1: The rotational frequency of Maclaurin spheroids as a function 
of their eccentricity (upper curve, the star indicates the point of onset of 
dynamical instability to a toroidal mode), and the Jacobi sequence (lower 
curve), after Chandrasekhar [1]. Also shown (thick dots) are numerical mod- 
els of quark stars calculated by us with the general relativistic spectral code 
of Gourgoulhon et al. [3]. The circle indicates that spheroid for which the 
marginally stable orbit grazes the equator (e = 0.834583). 
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Figure 2: The maximum frequency in stable circular orbits around the 
Maclaurin spheroids as a function of the spheroid eccentricity (continuous 
curves). The thin curve is the orbital frequency at the equator, eq. (2), 
the dashed portion corresponds to unstable orbits at the equator. The thick 
curve is the orbital frequency in the marginally stable orbit, eq. (3). The 
dotted portions of the curves correspond to spheroid eccentricities past the 
point of onset of dynamical instability, at e = 0.95289 [1]. Note that the 
curves are Newtonian, the speed of light does not enter eqs. (2), (3). Also 
shown (as thick dots) are the numerical models of quark stars presented in 
Fig. 1. 
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orbits around a quark star of (baryon) mass O.OIM©, modeled with the MIT- 
bag equation of state p = (p — po)c^/3. The actual values in Fig. 3 are for 
Po — 4.28 X lO^^g/cm^. The results for the Dey et al. equation of state 
[4] are similar, but the numerical values are higher as frequencies scale with 
the square root of density. The largest value of orbital frequency is obtained 
for the static model {e = 0, Q = 0), and is given by the Keplerian ex- 
pression a; (a) = ^AnGp/S ~ ^ifrGpo/S. For every e < the maximum 
stable orbital frequency is attained in a circular orbit on the equator, i.e., 
/ = (27r)~"^a;(a). For larger eccentricities 27r/ = <^(rms) = <^ms ^-nd the 
value of orbital frequency in the innermost stable orbit drops precipitously 
as e — > 1. The numerical results are seen to agree with the Newtonian ex- 
pressions of eqs. (2), (3). Clearly, general-relativistic effects are negligible 
for low-mass quark stars. 

Because in Fig. 3 the dependence on rotational frequency (27r)~^Q is 
shown, the curves are "double valued" (c.f. Fig. 1), for a given rotational 
frequency the larger value of orbital frequency corresponds to a spheroid/star 
of lower eccentricity, the smaller to the one for larger e. In this figure, the 
intersection of the curves for ujras and for a;(a) on the unstable branch (dotted 
curve) is a mirage caused by projection onto the Q axis, in fact, values of 
about 1250 Hz are attained for different values of e in the case of {27r)~^u}ms 
and of (27r)^^c<j(a). Compare Fig. 2, where it is clear that uj{a) and uj{rms) 
intersect only for e = Cg and e = 1. Dynamical instability to a toroidal 
deformation sets in at e = 0.95289 in the Newtonian theory [1], for this 
reason we do not extend the curves in Fig. 3 to low orbital frequencies. 

Finally, in Fig. 1 we plot the the rotational frequencies of the numerically 
computed stars, superimposed on Maclaurin's and Jacobi's sequence [1]. The 
circle corresponds to the appearance of the marginally stable orbit at e = Cg. 
Note that this occurs above the bifurcation point, at eccentricity ej = 0.8127, 
past which the Maclaurin spheroids are secularly unstable to deformation 
into a Jacobi ellipsoid [1]. However, this does not necessarily imply that 
the Newtonian marginally stable orbit, discovered here for the Maclaurin 
spheroid, is irrelevant to quark stars or other stars. Wc find that in our 
numerical models, the dependence of rotational frequency on eccentricity for 
higher mass stars (e.g., up to I.5M0) departs by no more than several percent 
from the curve of Fig. 1, at least up to its maximum. Similar departures are 
present in a post-Newtonian calculation of incompressible stellar models by 
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Figure 3: The maximum orbital frequency as a function of the rotational 
frequency of the Maclaurin spheroid. All symbols have the same meaning as 
in Fig. 2. 
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Shapiro and Zane [5], who find that the secular instabihty point progressively 
moves up to higher values of ej and to higher values of Q, until the bifurcation 
point reaches the maximum of the rotational frequency curve, ft{e). 

For a range of masses of compact stars, possibly existing in the real word, 
effects of general relativity may delay the onset of secular instability to non- 
axisymmetric deformations past the eccentricity where the essentially New- 
tonian marginally stable orbit (related to rotation-induced stellar flattening) 
appears. What is the minimum mass of a quark star for which ej > e^? Cal- 
culations in ref. [5] indicate that the bifurcation point occurs at ej — 0.835 
for GM/ {Rc^) = 0.013, where R is the radius of the non-rotating star of 
gravitational mass M. For the quark stars presented here, this would occur 
already for a star with M = O.O3M0, for which R = 3.2 km. It would appear 
then, that in the mass range ~ (O.OSMq, O.IMq) quark stars are well approx- 
imated by the Maclaurin spheroids and post-Newtonian effects may stabilize 
the quark star at eccentricities sufficiently high that the purely Newtonian 
marginally stable orbit is present outside the stellar surface, for e > e^. 

The expressions (2), (3), for the orbital frequencies around Maclaurin 
spheroids, and the presence of the marginally stable orbit [6] in their New- 
tonian theory are interesting in their own right. In addition, the gravity of 
Maclaurin spheroids provides a consistent Newtonian framework in which the 
infiuence of the marginally stable orbit on accretion fiow can be analytically 
investigated in two or three dimensions — up to now analytic work on accre- 
tion in black hole and neutron-star systems often relied on the convenient, 
but artificial, one dimensional "potentials" of the Paczyhski-Wiita type [7] . 

We note that the kHz QPO phenomenon (millisecond variability of X-ray 
fiux in bright Galactic sources [8]) has been interpreted as giving support to 
the presence of the general-relativistic marginally stable orbit in accreting 
neutron stars, and an application of the relativistic equivalent of Kepler's 
third law to the specific maximum values ~ 1 kHz of the observed frequen- 
cies in several such systems yields stellar masses ^ 2Mq in that interpre- 
tation [9]. But we saw that the same value of ujms/i'^'n') ~ 1 kHz can be 
attained in a (Newtonian) marginally stable orbit for a rapidly rotating star 
of rather low mass (e.g., lO'^M© or lO^^M©). Essentially, this is because 
in general relativity the gravitational radius = 2GM/c^ is a natural unit 
of length, and the characteristic frequency is ^JCM/Vg oc M~^, requiring a 
large mass for a low frequency, while for the Maclaurin spheroid the char- 
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acteristic unit of length is ae and the corresponding characteristic frequency 
is Uc = \JGM/ (ae)^. Since by mass conservation a?{l — e^)^/^ = const, 
ujI cx e~^(l — e^)^/^ for the Maclaurin spheroids and ujc vanishes for e = 1. 

We conclude that it would be premature to assign a mass value to a 
compact star based on observations of maximal orbital frequency alone, even 
if there is no doubt that the frequency is obtained in the marginally stable 
orbit. In a source where much of the luminosity is released in an accretion 
disk, it would be difficult in most cases to discriminate directly between a 
star of mass IAMq with a 10 km radius and, say, a O.IM© star of radius 
5 km. One indication of a low mass could be a lower photon flux during an 
X-ray burst, when the luminosity is thought to reach the Eddington limit 
value LEdd = 9 x 10^''(M/Mq) erg/s [10]. It could well be that many of the 
kHz QPO sources are in fact rapidly rotating low-mass quark stars. Fig. 3 
makes it clear that in such a case the maximum orbital frequency may have 
a value far lower than would be the case for a spherical star. 

Low-mass quark stars are well described by Maclaurin spheroids. This 
is the only known example of an accurate analytic description of a compact 
stellar remnant composed of matter at supranuclear density. We expect 
that in general, qualitative features of the exterior gravitational potential of 
Maclaurin spheroids may be relevant to a discussion of rapidly rotating stars. 
White dwarfs do not have uniform density, but if they rotate very rapidly, as 
they are expected to in cataclysmic variables, their orbital frequencies may 
turn out to be closer to those of the Maclaurin spheroids than to those of 
a point mass. Neutron stars are unstable below ~ O.IM0 while the results 
presented here are strictly applicable only to low-mass compact stars, but we 
speculate that the exterior metric of massive (M ~ Mq) neutron stars and 
quark stars in general relativity may share some features with the behavior 
uncovered here for Maclaurin spheroids. For example, in black holes frame 
dragging draws the marginally stable orbit in towards the axis of rotation, as 
the black hole spin increases, but it is known that for rapidly rotating massive 
quark stars rotation acts in the opposite sense: the marginally stable orbit is 
pushed out [11]. Hints of the same effect are observed in neutron star models 
at highest masses. We believe that the origin of this behavior is Newtonian, 
and point out that effects of stellar flattening, present already in Newtonian 
theory, counteract relativistic effects of frame dragging. 
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